Dynamics of particles and manifolds in a quenched random force field 
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O , We study the dynamics of a directed manifold of internal dimension _D in a d-dimensional ran- 

■ dom force field. We obtain an exact solution for d oo and a Hartree approximation for finite 

d. They yield a Flory-like roughness exponent C, and a non trivial anomalous difi^usion exponent 
V continuously dependent on the ratio Qt / of divergence- free {qt) to potential (pl) disorder 
strength. For the particle {D = 0) our results agree with previous order RG calculations. The 
^ I time-translational invariant dynamics for > smoothly crosses over to the previously studied 

C . ultrametric aging solution in the potential case. 
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The problem of the diffusion of a particle in a d-dimensional static random force field has been often addressed in 
the past decade It is of interest in problems like tracer diffusion in porous media |^,^ or turbulent flows 

H ■ Typically the force is the sum of a potential random force (longitudinal in fc-space) of strength gi^ and a divergence-free 
\ (transverse) one of strength g^. The diffusion is anomalous [0-|6j (r ~ with v ^ 1/2) if either (i) the dimension d 
P3 ' of space is low (d < 2), or (ii) the force field has long-range correlations {F{Q)F{r) ~ with a < 2). While this 
O is true for weak gaussian disorder, other anomalous behaviors are possible for disorder distributions with tails 

Apart from particular cases (such as the long-range correlated potential case [Q, including d = 1 [^), most of the 
above results were obtained by perturbation theory in the disorder around the thermal value i^th = 1/2, which holds 
in the quadrant d > 2 and a > 2. An expansion in powers of e' = 2 — d was performed for short range (SR) disorder 
[01, and one in powers of e = 2 — a for long range (LR) correlated disorder |^,^ (and a double e, e' expansion ||^ near 
the point d — 2, a — 2). Superdiffusive behavior was obtained for gi/gi^ 3> 1 and subdiffusive slow dynamics for 
^ f/x/ffL ^ 1- A line of fixed points was obtained up to order e^, with a continuously variable exponent in the case 
Cn| ' of LR correlations (d > 2, a < 2), parametrized by the ratio ^x/ffL- The divergence-free case = is simpler, since 
the value of v is equal to the Flory value obtained by purely dimensional arguments, as discussed below. 

On the other hand, out of equilibrium relaxational (potential) dynamics has recently been studied extensively, in 
the context of spin glass models |^ as well as particles and manifolds in random potentials. Analytical 

"j^ solutions were obtained |^ for the long-time dynamics in mean-field models, showing explicitly their glassy features 
such as slow (aging) dynamics, and the breakdown of time-translation invariance (TTI) and of fluctuation-dissipation 
I ■ theorems (FDT) . Whether these glassy properties of potential dynamics survive to a small non-potential perturbation 
■ O [ was recently addressed and it was found that, while the dynamics becomes TTI in most cases, some glassy 
^ ■ features survive. Similarly, a glassy state with stationary dissipation was proposed for driven vortex systems. 

In this paper we thus reconsider the dynamics of a particle [D = 0) in a d-dimensional random force field as well 
as its extension to a directed elastic manifold of internal dimension D. We obtain the diffusion exponent r and 
(for Z? > 0) the roughness exponent r ^ x'' where x is the internal coordinate. Our results are exact for d — > oo. 
^ . This mean-field limit is non trivial if transverse and longitudinal parts are scaled properly. It extends the previous 
?H results obtained via RG since it allows to reach nonperturbatively arbitrary values of the ratio gT / gh (that 
the RG could not reach because of runaway flows for = 0). We can now explore systematically the interesting 
potential limit 5t ^ and we prove the remarkable property of continuity of the mean-field dynamic solution when 
aging is approached for ^ 0. This property was proposed in |]l3| and verified in numerical solutions of the p-spin 
model. We suggest to take systematically advantage of it, in order to improve numerical simulations by adding a 
small nonpotential part, waiting till TTI is attained, and taking the nonpotential part to zero after the large size 
limit. We find that the finite-d Hartree approximation compares well with the RG. 

We study the Langevin dynamics of a directed elastic manifold, of internal dimension D, parametrized by a d- 
component field r{x): 

dtVaix^t) = Vlra{x,t) + Fa[r{x,t),x] +ria{x,t) , (1) 
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where rj is the gaussian thermal noise, satisfying {r]a{x,t)i]i3{x' ,t')) = 2TSai3S{t — t')6^{x — x'), and Fa[r,x] is the 
gaussian quenched random force with a correlator Fa[r,x]Fp[r',x'] — Aa^(r — r')6^{x — x') which depends only on 
the distance |r — r'|. We denote by angular brackets averages over the thermal noise and by an overbar averages over 
the random force. The force has both a potential (L) and a divergence-free (T) part, which contribute separately to 
the correlator: 

A„^(r) - -dr^dr,^i.{r) - [Sa^p ^r, - dr^dr,"^ ^T{r) . (2) 



We study the case in which the correlations at large scale behave as power laws: the Fourier transforms of the 



correlators with respect to r behave, for small values of their argument K, as Aap{K) ~ K"^ ''■{gi^P^i^{K)+gTPaBi^))j 



where P^p{K) = KaKp/K"^ is the longitudinal projector and P'^p{K) = 1 — P^p(K) is the transverse one. We 
compute the mean-squared displacement Dxx'it,t') and the correlation Bxx'{t,t'), respectively defined by I?^;^:' (i, t') = 
^((r(x,i) — r{x' ,t')Y) and Bxx'{t,t') = i((r(x,i) —r{x,t')) ■ {r{x',t) — r{x' ,t'))) . The mean-square displacement at 
equal times Dxx'{t,t) ^ \x — x'p'" measures the roughness of the manifold, while the correlation at equal space 
B{t,t') — Bxx{t,t') measures the diffusion of a single monomer. When TTI holds, the diffusion exponent v is defined 
by B(t, t') ~ |< - t']^". In the absence of disorder one has C = Cth = (2 - D)/2, v = i^th = (2 - D)IA for D <2, and 
bounded motion and fluctuations for D > 2. Before proceeding, let us consider the problem by Flory-like dimensional 
arguments. If all terms in (|l]) scaled the same way one would have r/t ^ rjx^ ^ r^^x^^/"^ with 7 — a/2 for d > a. 
This yields 

Cf - - t4-^ • (3) 

4(1 + 7) 

It turns out that, for d — > 00, C = Cf is exact but ly — holds only in the transverse case. 

One can take the large-d limit directly in the De Dominicis-Janssen functional, performing a gaussian decoupling 
of r{x, t) and the conjugate field t). We define: 

AtW = -V^ir^d) , AlW - -d Vi^ir^d) . (4) 

One can also study in any d the self-consistent Hartree approximation where r{x,t) and f{x,t) are treated as 
Gaussian variables. It amounts to replace everywhere, in the equations and solutions obtained below for d — 00, 
14,T by Vl,t defined by V{z) = J d'^(?i/(27rz)''/2 exp(-0V(2z))y(0Vd) (see, e.g., appendix A of 0). The coupled 
dynamical equations also involve the response Rxx,{t,t') — l/dS{r°'{x,t)) /6f/^{x' ,t')\f=o where f{x',t') is a small 
perturbation applied at point x' at time t'. The response and correlations are isotropic, {t, t') — {5ai3 / d)Rxx' {t, t'). 
We choose a space-homogeneous initial condition Bxx'{0,0) = Bx-x'{0,0) and denote by Bk and Rk the Fourier 
transforms with respect to a: — x', and by tilde equal space {x — x') two-time functions (e.g., B = J^, Bk). We use 
the shorthand Jf. = J d^k/i2n)^ and = J dui/{2TT). With Rk{t,t) = and Rk{t,t~) — 1, the dynamical equations 
read 

.2) 



dtRk{t,t') = -k^Rk{t,t')- / dsJ:{t,s){Rk{t,t')-Rk{s,t')) + 6{t-t') , (5) 

dtBk{t,t') = -k^ Bk{t,t') + 2T+ f ds D{t,s) Rk{t,s)- f ds D{t, s) Rk{t' , s) 

Jo JO 

t 

dsi:{t,s) {Bk{t,s)-Bk{s,t') + Bk(t,t'))-4TRk{t',t) , (6) 







where D{t, s) = W{iB{t, s)) and S(i, s) = -W^'{B{t, s))R{t, s). We define Vi = Vl + [l - l/d)VT and V2 = Vt, and 
consider Viib) = gib^-^/{2{l - 7)), ^2(6) = g2^''"^/(2(l - 7)) where gi = gi^ + (1 - l/d)gT, 92 = 5l, and 7 = a/2 
for large d, thus generahzing the potential case studied in pl| , p^ . We define gL,T — —.9L,TCd, 7/(2(1 — 7)), where 
C(i,7 = / (r^/d)"'^"''' exp(in • r) with — 1. As shown below, the case of LR correlations corresponds (for d 00) to 
7 > 2/(2^^- D) for D <2, and any 7 for D > 2. 

When (7t 7^ one can assume that a TTI solution is reached at large t' , Bk{t, t') = Bk{t—t') and Rk{t, t') ~ Rk{t~t'). 
We checked this property numerically for £> = 0. However we stress that it is not certain that it should hold for any 
nonpotential dynamics p6|. Using TTI one can recast the above equations in the form 
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Rk{(^) = — ,o , ^,n^ ' -^fci'^) = T? 27r(5(w) - 



/+OC /• + 00 

dtexp{iujt)V2ibit))r{t), D{uj) = A j dt exp{iLut)Vl{b{t)) , (7) 



where L'fc = J^,{D{uj')Rk{-uj') - Y.{uj')Bk{uj')). The correlation Cfc(a') = l/d{{r{k,uj).r{~k, -uj)) thus obeys 
2Cfe(ci;) = (4T + - icj + fc^ + S(0) - S(a;)p. The displacement can then be computed using Dx{t) = 

2 /j, — cos(fca; + ujt))Ck{i^)- In the large time limit we propose 

R{t) Rt-f" , B{t) Bt"^" , (8) 

and define the FDT violation factor X[B] by R{t) = X[B{t)]dtB{t). Equation (^) implies the small frequency 
behaviors R{uj) - RT[1 - ^j]{-iuj)''-^ for ^ < 1 or R{uj) - R{0) - RT[1 - ^j]{-iuj)''-'^ for fi > 1 and B{lu) - 2BT[1 + 

2iy] Reiiu)-^-^" . It also implies D{uj) ~ 8giB'^T[l - 2u-f] Reiiuf"^-^ and 2(0) - S](w) 4,g2iRB~^^+^^T[l - 

— 2z/(7 + l)]{—iu;)'^~-^^^'^^'^^^\ One must have < ^ — 1 + 2^(7 + 1) < 1 in order to be able to neglect the —iuj in 
the Green functions and to have a finite S(0). We shall find out that this is always the case for 52 > 0. The equation 
for R = /j, Rk yields, at small ut and for D < 2, 

i?M -v4z3(E(0)-SM)^ , (9) 

with Ad 7r/(2^7r^/2p[£)/2] sin(7rD/2)). For D > 2, R{Q) is finite and one has instead R{uj) - R{0) ~ Ad{T,{0) - 
T,{lu))^~ . This gives the relation between exponents: 

l-Ai = 2K7 + l)|^ • (10) 

This value for fi is consistent with neglecting the —ilu term in (|^), provided that < < j/t, where i^x = (4 — 
i:')/(4(l + 7)). It implies X[B] - XB'" with with X = R/ {2v)B-^^+'^'^^ and w = -1 + (1 + 7)|5§. It also yields 
the following relation between amplitudes: 

,0 . (-4g27r[2z.(7 + l)4^])'^ ,,,, 

(2l/X) 2 ^ An TT-rT- ■ (11) 

^ ' ° r[2K7 + i)|5g] ^ ^ 

Similarly the equation for B — Bk yields 

Bic.) ~ AnDic.)\m - S(a;)|^ '^^!'^X^ ' (^2) 

sm(0j 

where </) = — f (a* ^ 1 + 2!/(7 + 1)) is the phase of S(0) — S(cj). One can check that the exponent relation ( p^ is now 
automatically satisfied. The resulting relation between amplitudes, together with the relation (pT|), yields our main 
result, i.e., the equation which determines the exponent v. 

9l=h( n )- ^[^^ + 1)1^] r[l-27Hsin(^7^)sin(f2Kl+7)|^) 

51 |r[2K7 + l)45%]l r[l + 2H7sin(^^) sin(f2Kl+7)4^) ' ^ ' 

where, we recall, 52/31 = .9l/(5l + .9t) when d — > cxd. The function h is a monotonic function of u, such that 
h{iy — 0,15,7) — 1 and /i(i^ = ^^,0,^) — where j^f is given by (H). Thus we find that in the transverse case 
ly — I'T = I'F and that ly spans the interval Q < v < i^t when the ratio 52/<?i varies from 1 (potential case) to 
(divergence free case). These results hold in the LR case 7 < 7c. The limit 7 ^ 7^7 (for which /i(i^,D, 7c) = 1) is 
studied below. One also easily obtains the roughness exponent C- From Cfc(w) and Eq. (|^) one finds that i ~ a;^ with 
z = 2v/v^ and r ^ . Thus C, ^ zv = 2i^T — Cf- This is the Flory value, as announced, and interestingly, is found 
not to depend on the precise nature of the random force (though the physics is quite different as the ratio gh/gT 
varies). For finite d our approach also yields a Hartrec approximation which amounts to: (i) replace in ( p^ gi by 
.gL + (1 ^ i)9t ^^'^ 52 by ^l; (h) setting 7 = min(a, d)/2 (7 = d/2 corresponds to the Hartree SR case and 7 < d/2 
the Hartree LR case since the model studied here and in mM is a random force model. 
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Let us compare our result for v (exact when d — > oo) with the RG calculations for the particle [D = 0) with 
7 = 1 — e in any d. The RG calculation was carried out upto order e in ||^ and upto order in |3|. In our notation, 
it reads 

1 d-l-dk ^ ,4 + 2d-5d^ + d^ + dk(4:-2d + d^) , , 

Y. = ^- -W^' + ' * 4i2-d)id-iyd ' 



where k = gL/gr, while, from (^^ and using the above Hartree values one finds 

1 _ . .91 - 2g2 g2.gie^ _ d-l-dk e^dk{l - d - dk) 

2v~ 2{g,-g,f A{gi - g2? ^ 2(d-l)'+ A(d - If ' ^ ' 

For d — > oo, 1/ from (|l3|) and (|l|) coincide. V is determined by the bare ratio gx/gL, which indeed is unrcnormalized 
in the RG. Furthermore the finite d Hartree approximation is found to be exact to first order in e = 1 — 7 for all 
d > 2. Thus it is a good approximation in the Hartree LR case and is exact for d = 1 since it yields ly — 0. Though 
it does predict correctly the special point c? = 2, 7 = 1 it is not able to reproduce correctly the delicate crossover |^ 
around this point. The RG approach of can be generalized to the manifold, and for < D < 2 one can similarly 
expand around jdD) = 2/(2 - D) (LR) (and around 4(D) = 4/(2 - D) (SR)). The resuh (|l|) leads to the following 
formula, exact for d —> 00, with 7 = 2/(2 — D) — e to first order in e: 

2-D (2 ^ D)^ (4 - D)g2 - 2g, 

It would be interesting to check whether the Hartree approximation is also exact for any D < 2, d to 0(e). 

Up to now we have studied the limit 7 lc{D) keeping the ratio gi^/gT > fixed, in which case v vth- In 
fact the limit 7 — > 7c is very non uniform, and therefore much richer. For D < 2 one finds that any value of v with 
< < Uth can be reached if one takes the limit 7 —> lc{D) with the ratio g^/egij fixed. For the particle, this is 
perfectly consistent with the RG result obtained in any d, which in our notations read 1/v = I + 5L + 0{gi^) 
(note that it has a nice d ^ 00 limit and that there is a conjecture that all higher order terms vanish). For 
2 < D < 4 one finds that the limit 7 ^ 00 is non trivial with a function of gT/gh- Finally the limit D = 4 — e is 
also very non uniform and any value of v can be reached if gx/e is kept fixed as e ^ 0. 

We now study the potential limit gx —^ and show on the case D = that one goes continuously to the aging 
dynamics which holds for g^ — 0. This is best illustrated using the formalism of Ref. |l^. Both in the case of a 
solution satisfying (gx > 0) and in the asymptotic aging solution {gj^ = 0), one defines the FDT violation factor X[B] 
and the triangular relation f{B,B') by 

X[B]= lim R{t,t')/dt'B{t,t') , J{B',B)= lim B{s,t') . (17) 



t>t'^oo t>s>t' 

Blt,t') = B>B(t,s) = B' 



The dynamical equations are naturally written in terms of the linear susceptibility x(t, t') = J*, dsR{t, s), which takes 

the form xiB) = Jq dB' X{B'), and of the "anomaly" M[B] = 4 /^f dS' V"{B')X{B'). The equation for R only 
contains V2 and is thus identical to the case gx = and reads, in the long time limit, like Eq. (5.3) of Ref. pl| : 

« ^ = -1 - X[B] dB' A\q'{B')X[B'] + Jj dB'Wi'{B')X[B'] {x[l{B' , B)] - x[B]) . (18) 

We used B{t, 0) — > 00. This equation is coupled to a similar (more cumbersome) equation for B that also involves Vi. 
We specialize to the long range case 7 < 1. The integrals are naively divergent at large B since lims^oo x[B] — 00 
In the equation for B similar divergences occur. However, at large B, one finds x[B] ^ 2X/(1 + 7)^^^+''')/^, as 
well as f{B', B) = {B^ — B'^ Y'^ . The divergent parts cancel and one recovers the equation for the amplitudes (|ll| ) 
(independent of small B details). This solution interestingly shows that when — > as t/r — *■ the triangular relation 
continuously becomes f{B',B) — max(_B',_B), i.e. the known ultrametric solution for the potential case. Similarly 
X[B] ~ XB^'^-i)/^ tends continuously to the resuh of 

In the transverse case gi^ = for the particle D — Q our method yields the Flory estimate v — v-p = 1/(1 + 
min(a,d)/2), which turns out to be exact (for v-p > 1/2). This can either be seen perturbatively to all orders in the 
field theory j^, or, more physically, using that the fact that the stationary measure is exactly given by Pgtat — const. 
The Derrida-Luck method 0,|l^ then produces the result (and also yields rigorously the velocity-force law v ^ F'^ 
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with (j> = 1) as we now show. One introduces a medium with random forces, but periodic with a large period A. An 
exact formula for the velocity in this medium is wa = J dxF{x)Pststt where Pstat(a;) is a stationary measure (satistying 
= V^Pstat — di{FiPstat)) with periodic boundary conditions. Using that Pgtat = l/A'^ one finds that the mean- 
squared velocity v\ ^ \-min{a,d) -^^ ^^lat periodic medium (equivalently finite torus of volume A''). This exact scaling 
relation ensures that v — vp. For the manifold D > 0, however, the stationary measure is not known (the naive 
extension Pstat ~ exp(— / da;(Vr)^/2r) is not stationary). Thus the Flory values obtained here for the manifold in a 
divergence free flow may very well be exact only for d — > oo. 

One can also extend the arguments of Ref. [|| to the diffusion of an oriented (e.g. by an external field) manifold in 
a divergence-free turbulent flow. If the random force, proportional to the local velocity, can be considered as static 
on some time scales, with LR gaussian correlations v{0)v{r) ^ ^^/s+t^/s arises from the intermittency corrections 
to Kolmogorov scaling), the result v = vt for the corresponding diffusion exponent yields i/ = 3(4 — £))/(2(4 — /i)), 
which can be considered as a generalization of the Richardson's law (corrected by intermittency). 

In conclusion, we have studied the problem of the dynamics of a manifold or a particle in a gaussian random force 
flow. We have obtained the roughness exponent C and v exactly for d — >• cxd. Our study establishes a bridge between 
several problems (aging and nonpotential dynamics) and techniques (large d and RG approaches). 
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